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ABSTRACT 

Aims. We seek to provide additional tests of the line formation of theoretical 3D solar photosphere models. In particular, we set out 
to test the spatially-resolved line formation at several viewing angles, from the solar disk-centre to the limb and focusing on atomic 
oxygen lines. The purpose of these tests is to provide additional information on whether the 3D model is suitable to derive the solar 
oxygen abundance. We also aim to empirically constrain the NUTE recipes for neutral hydrogen collisions, using the spatially-resolved 
observations of the O i 777 nm lines. 

Methods. Using the Swedish 1-m Solar Telescope we obtained high- spatial-resolution observations of five atomic oxygen lines (as 
well as several lines for other species, mainly Fe i) for five positions on the solar disk. These observations have a high spatial (sub- 
arcsecond) and spectral resolution, and a continuum intensity contrast up to 9% at 615 nm. The theoretical line profiles were computed 
using the 3D model, with a full 3D NLTE treatment for oxygen and LTE for the other lines. 

Results. At disk-centre we find an excellent agreement between predicted and observed line shifts, strengths, FWHM and asymme- 
tries. At other viewing angles the agreement is also good, but the smaller continuum intensity contrast makes a quantitative comparison 
harder. We use the disk-centre observations we constrain ^Sh, the scaling factor for the efficiency of collisions with neutral hydrogen. 
We find that 5h = 1 provides the best match to the observations, although this method is not as robust as the centre-to-limb line 
variations to constrain ^Sh- 

Conclusions. Overall there is a very good agreement between predicted and observed line properties over the solar granulation. This 
further reinforces the view that the 3D model is realistic and a reliable tool to derive the solar oxygen abundance. 

Key words. Sun: granulation - Line: formation - Sun: photosphere - Techniques: spectroscopic - Techniques: high angular resolution 



1. Introduction 

One of the significant problems in modeling photospheres is 
how to accurately treat convection. Three-dimensional time- 
dependent hydrodynamical simulations of convection in solar 
and stellar photospheres are a major leap forward in the treat- 
ment of convection and realistic physics. They challenge the 
classical approximations: plane-parallel photospheres, static hy- 
drostatic equilibrium and the mixing-length theory of convec- 
tion. 



Employing the family of convection simulations of Stein & 

frdlund (1998), the series of papers started by |Asplund et al. 
00) has undertaken a systematic study of solar line formation 
and elemental abundances. The solar abundances derived from 



these models ( Asplund et al. 2005a 2009) represent a down- 
ward revision of the solar metalicity by almost a factor of two. 
The significance of this revision echoes over many fields of as- 
trophysics, because the Sun is a reference in abundance studies. 
In solar physics it has sparked much debate, mostly because the 
excellent agreement between the solar interior models and mea- 
surements from helioseismology is lost when using the revised 
solar abundance s ([Bahcall et al.||2005| ). Despite numerous tries 
(see review of |Basu & Antia||2008| ), no work has successfully 



reconciled this excellent agreement with the new solar abun- 
dances. 

Realistic physics, absence of free-parameters governing tur- 
bulence and the excellent observational agreement of predicted 
spectral line shapes ( [Asplund et al.|2000) are some of the strong 
points in favour of the 3D models. Yet the disagreement with so- 
lar interior models and helioseismology has driven the question: 
are the 3D models suitable to abundance analysis? 

This work is part of an effort to addresses several observa- 
tional tests in order to answer this question. In particular, we 
use high- spatial-resolution observations to test the 3D model 
and line formation over the solar granulation, and at diflTerent 
viewing angles. A robust test of the models can be established 
by comparing the observed and predicted properties of spectral 
lines over the range of temperatures and velocities found in the 
solar granulation. Looking at diflferent angles (or positions on 
the solar disk) one can also probe higher photospheric layers and 
non- vertical velocity components. Such testing has been briefly 
discussed by [Asplund] ( [2005[ ). The series of papers started by 
Asplund et al. ( [2000 ) shows that the 3D models reproduce the 
observed mean spatial line profiles very well. In this work we 
aim at resolving those mean line profiles and test the 3D models 
for a few lines. 
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We take a special interest in oxygen lines. Oxygen is piv- 
otal in the solar abundance revision. Not only because it is 
abundant but also because abundances of several other elements 
are determined by t heir ratios to oxygen. After the work of 
[Asplund et al. ( [2004| ) several new determinations of the solar 
photospheric abundance h ave been made: either using a differ- 
ent family of 3D models (jCaffau et al. ^'2Q08^, "^res^SOOS ) or 
different methods ( Socas-Navarro & No rton 2007; Centeno & 
Socas-Navarro|2008| ). The discussion provided by these studies 



prompts us to revisit the spectral indicators of oxygen. We pro- 
pose to study some of the oxygen lines used in the abundance 
analysis. High- spatial-resolution observations of five atomic 
oxygen lines were obtained for this purpose. Observations were 
also obtained for other lines close to the target oxygen lines, and 
we have included some of them in our analysis (mainly Fe i lines, 

see Table |2]). 

This work is closely related to Pereira et"aLl ([2009 , hereafter 



Paper which focuses on the centre-to-limb variation of oxy- 



gen lines. We make use of the same observations, but here we 
focus on high-spatial-resolution, whereas "Paper 11' studies the 
spatially and temporally averaged spectra as a function of view- 
ing angle. The discussion here is primarily focused on the disk- 
centre granulation, because at disk-centre the contrast in inten- 
sity between granules and intergranular regions is highest. This 
greater contrast gives more information about the variation of 
spectral lines than at other viewing angles. 

To study the line parameters we employ a similar methodol- 
ogy as |Kiselmanl ( |199 4 ). We study the variation of line strengths, 
shapes and shifts over the observed granulation and compare it 
with the 3D model. For the O i 777 nm triplet lines, known to 
suffer departures from LTE (g.g. [Eriksson & Toftj 1979 [[Asplund] 
2005[[Fabbian et al.[2009[ ), we seek to find whether the spatially- 



resolved disk-centre observations can constrain the NLTE line 
formation, as has previously been done using the centre-to-limb 
variations ( Allende Prieto et al.|2004 Paper II j ). 

The outline of this paper is as follows. In Sect. [2| we detail 
the observations. In Sect. [3] we describe the data reduction proce- 
dure. Next, in Sect. [4] we provide details of the 3D model used, 
the line formation, atomic data and how the model results were 
degraded to the observed conditions. In Sect. [5] we present the 
results and discuss their significance. Finally, we establish some 
conclusions in Sect. [6] 

2. Observations 

2.1. Telescope, instruments and programme 

The observations were acquired using the Swedish 1-m Solar 
Telescope (SST) ( [S charmer et al.[[2003a[ ) on Roque de Los 
Muchachos, La Palma. The observations took place from the 
14th till the 29th of May, 2007. 

Our aim was to obtain high-quality spatially -resolved spectra 
across the solar disk: from the disk-centre to the solar limb, with 
a high spatial and spectral resolution. We made use of TRIPPEL: 
the TRI-Port Polarimetric Echelle-Littrow spectrograpl|^ The 
slit width is 25 yum, which translates to 0.11 arcsec. The spec- 
tral resolution (A/AA) is ^200000. TRIPPEL allows the simul- 
taneous observation of three wavelength regions (spectrograph 
ports). To obtain spatially-resolved spectra along the slit, there 
is no predisperser. All spectral orders are overlapped and order 
selection is done by using filters. 

The spectrograph cameras used were two MegaPlus 1.6 cam- 
eras (615 and 777 nm windows) and one MegaPlus II esl603 

^ See http://dubshen.astro.su.se/wiki/index.php/TRIPPEL_spectrograph. 



camera (630 nm window). In addition to the spectrograph and 
its three cameras, we operated four imaging cameras: three slit- 
jaw image cameras (at the wavelengths of 630, 694 and 773 nm, 
synchronized with the spectrograph's cameras) and one 'blue' 
camera with a Can H filter, using a dichroic beam splitter for 
A < 500 nm. The slit-jaw cameras allowed for structure iden- 
tification; the Ca H camera was used as a monitor of magnetic 
activity. The exposure times used were in the range of 30 - 265 
ms, depending on the filters and the position on the solar disk. 
With neutral density filters to ensure identical exposure times, 
the slit-jaw cameras were synchronized with the spectrograph's 
cameras. 

Our main target being oxygen lines, the filters for the three 
wavelength regions of the spectrograph were the following: 
615.40, 630.08 and 777.52 nm (filter transmission windows 
are between 5.1-6.8 nm, centered on the filter wavelength). 
Hereafter referred as 615, 630 and 777 regions, they correspond 
respectively to the spectrograph orders 37, 36 and 29. These 
wavelength regions allowed us to observe the 1 615.81 nm line, 
the [Oi] 630.03 nm line and the three Oi triplet lines around 
777 nm. A summary of the lines of studied in these three regions 
is given in Table [2] 

High-spatial and spectral resolution observations of these 
lines were obtained for several positions on the solar disk, from 
the centre to the limb. The positions in the disk are defined by 
yu = cos 6, where is the heliocentric viewing angle. Because 
the 3D models used do not include magnetic fields, our aim in 
this work is to study the quiet Sun. While observing, efforts have 
been made to avoid regions with obvious magnetic fields (e.g. by 
using the Ca H camera). However, magnetic fields in the Sun can 
never be completely avoided. 

In total more than 890 000 images were taken. The data were 
taken in a variety of seeing conditions, hence the same quality 
was not possible for all observing runs. Extensive use was made 
of the SST's Adaptive Optics system (Scharmer et al."2003b|). 
Some difficulties of acquiring AO lock were experienced close 
to the limb where the contrast of solar features is usually low. 
In some cases (notably jj, = 0.4) targets with higher contrast, 
like facular points, were chosen to facilitate locking. These struc- 
tures were presumably magnetic features and the procedure thus 
somewhat compromised the effort to observe as quiet sun as 
possible. Another systematic difference between datasets results 
from the fact that the disk-centre observations had higher pri- 
ority. Thus more time was spent on them and they happened to 
catch the moments of best seeing. 

Details about the data and the coordinates in the Stonyhurst 
system ( Cortie|18"97] ) can be found in Table [Tl and the observed 
locations on the solar disk are shown in Fig^]T] For the 777 nm 
camera at ^ 0.8 the dataset from the 25th of May was used. 
This was due to problems with a misaligned 777 nm filter for the 
27th May set (the adopted set for the other cameras). 

2.2. Pointing accuracy 

To minimise the spread in ji for each data set, a position on the 
disk was chosen so that the slit was parallel to the closest part of 
the limb. The field of view is rotating in the laboratory frame and 
thus the position angle of the selected locations changed during 
the day. This explains the scattered appearance of the points in 
Fig.[l] 

The SST's blind pointing is not precise enough for yu to be de- 
termined with the required precision close to the limb. We took 
advantage of the Dutch Open Telescope (DOT) to help with the 
coordinates. With the kind assistance of P. Siitterlin, the DOT 
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Fig. 1. Position of the observations on the solar disk, in helio- 
graphic Stonyhurst coordinates. The figure has been centered in 
the heliographic latitude of the observer -1.8° for our dates). 



Table 1. Observing runs selected for analysis. Coordinates are 
given in the Heliographic Stonyhurst system. The spread in [i is 
given by the approximate spatial coverage of the slit. The A/rms 
is given as the maximum continuum contrast of a single frame 
in the 615 nm wavelength region. 



UT Date 








Coordinates 


A/rms [%] 


20/5/2007 08:00 


0.999 


+ 


0.001 


0.2S, 0.5E 


9.2 


25/5/2007 08:30 


0.816 


+ 


0.010 


26.7N, 22. IE 


5.4 


27/5/2007 15:30 


0.793 


+ 


0.012 


1.3N, 37.5W 


6.0 


25/5/2007 09:56 


0.608 


+ 


0.020 


51.1S,3.7E 


4.7 


20/5/2007 09:27 


0.424 


+ 


0.024 


29.9S,61.4W 


5.2 


25/5/2007 07:40 


0.197 


+ 


0.003 


18.0S, 82.4E 


3.4 



was pointed at the same solar region as the SST by compar- 
ing the Can H filtergram view of the two telescopes. The he- 
liographic coordinates were thus read from the DOT. The spread 
in [I in Table [T] is due to the range of [i covered by the slit. 

Note that the outermost images around = 0.2 have a very 
well-defined ji because the limb was in the field of view of the 
slit-jaw images, which allowed precise measurements. This ex- 
plains the small spread in yu for that position. 



3. Data reduction 

3.1. Dark currents, flat fields and spectrograph distortions 

The first step in the reduction was the subtraction of the dark 
currents and division by a normalized flat field image. Darks 
were obtained by blocking the sunlight. Because of the features 
on the solar surface and the presence of the spectral lines, the 
task of obtaining a uniformly exposed flat field needs a more 
elaborate approach. The telescope is operated in a special mode 
where it scans randomly quiet areas at the disk-centre. A ref- 
erence disk-centre spectrogram is obtained by averaging these 
images. The spectral lines are removed by dividing this spatially 
uniform spectrogram by its mean spectrum (spatial average), and 



a flat field image is obtainecj^ This disk-centre mean spectrum 
is also used for wavelength, intensity and scattered light calibra- 
tions, as detailed below. 

Once corrected for dark currents and flat fielded, the images 
need to be corrected for the spectrograph's distortions: smile and 
keystone. The aim is to rectify the spectrogram, to have the dis- 
persion exactly perpendicular to the pixel rows. Smile curves the 
spectral lines across the spatial axis, whereas keystone curves 
them across the spectral axis. The distortion is computed by fit- 
ting a fourth order polynomial to the line centres of a few lines 
along the slit (in a spatially uniform image). This information is 
used to compute the transformation to invert the distortion and 
interpolate the spectrogram into a rectified version. The spectro- 
gram is also slightly cropped to remove dark borders. 

3.2. Wavelengtti and intensity calibrations 

Wavelength calibration is made by comparison with the ref- 
erence FTS atlas of |Brault & Neckel] ([1987 ). Using the disk- 
centre mean spectrum obtained from the flat fielding procedure, 
for each wavelength region a few lines are identified in the 
mean spectrum and in the FTS atlas. Their line centres are com- 
puted and a polynomial is fitted to the dispersion relation (wave- 
length/pixel number) - yielding a wavelength scale for our im- 
ages. The precision of this wavelength calibration has proven to 
be more than enough for our purposes. When computing the line 
shifts across the solar granulation, velocity measurements are 
made relative to the reference mean spectrum (the spatial and 
temporal mean of the images at a given ji) - thus a very pre- 
cise absolute wavelength calibration is not required. Using the 
telluric lines in the 630 nm region to calibrate the wavelengths 
more robustly we find a very good agreement (< 40 m s"^) with 
our 'normal' procedure. 

All spectral intensity measurements are made relative to the 
observed disk-centre intensity. In the flat fielding method de- 
scribed above, a spatially and temporally averaged reference 
disk-centre spectrum (hereafter referred as 'mean spectrum') is 
obtained. The intensity of each spectrogram observed at a given 
position on the solar disk is then divided by the (continuum) in- 
tensity of the mean spectrum. 

For the same exposure times, the exposure level varies with 
time-dependent factors, most importantly airmass. Because the 
flat fields are also used to calibrate the intensity, it is important to 
have the exposure levels of flat fields and observations as close 
as possible. To minimize for the time-dependent factors the ob- 
servations and the mean spectrum are obtained in close inter- 
vals (20-30 min). These intervals are usually enough to keep the 
same exposure level between observations and flat-fields (for the 
mean spectrum). However, in the early morning and late after- 
noon when the airmass varies quickly with time it is necessary 
to correct for the diff'erent airmasses. Our approach is empirical: 
using several thousand slit-jaw images for the given airmass in- 
terval, we fit the terrestrial atmospheric absorption coefficient r© 
from the logarithmic intensity v^-. airmass relation. 

We assume the terrestrial atmospheric extinction to be given 

by: 

/ = /o exp(-aTe), (1) 

where a is the airmass. The airmass correction factor is then the 
ratio of flat fields/observations intensity: 

/factor = exp(-T©[aff - Gi]), (2) 



^ For the spectrograph's cameras. For the imaging cameras the flat- 
field is simply the normalized averaged of the many disk-centre images. 
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where is the mean airmass when the flat fields were taken 
and Gi the airmass when each image was taken. The intensities 
from slit-jaw images are then used to derive r© empirically. With 
a linear fit to the log(/) vs. airmass relation we obtain r© and 
log(/o). Due to the absence of a dedicated slit-jaw camera, for 
the 615 nm region r© was extrapolated in wavelength from the 
other regions. 

In some circumstances the camera shutters introduce a 
slightly diff'erent exposure between odd and even frames. This 
eff'ect is corrected by scaling the exposure level of each odd/even 
frame by the ratio of the mean intensity to odd/even intensity. 

3.3. Correction for stray light in the spectrograph 

Light scattered off' the grating of the TRIPPEL spectrograph as 
well as from other internal sources introduces straylight in the 
spectra. In our treatment for straylight, we have assumed that it 
is constant along the slit and is proportional to the mean inten- 
sity of light entering the spectrograph. Our treatment was based 
on routines used in AUende Prieto et al^ ( |2QQ4 ), that were kindly 
provided by C. Allende Prieto. The relative level of straylight, 
given as a fraction of the continuum intensity, is computed us- 
ing the FTS atlas and the calibration disk-centre mean spectrum. 
First the mean spectrum is normalized. Assuming that the spec- 
trograph's instrumental profile is Gaussian, we derive the stray- 
light level and the resolution relative to the FTS atlas by compar- 
ing our mean spectrum with the atlas (convolved with a Gaussian 
for an arbitrary resolution). The approach searches for the best 
fitting values of the straylight and the FWHM of the convolv- 
ing Gaussian. Using an iterative minimization, the best fit- 
ting pair is found. This straylight and resolution represent the 
parameters that best 'map' the FTS atlas into our mean spec- 
trum. This procedure is repeated for each wavelength region and 
each data subset. We find a straylight level of ^6%, which is 
in agreement wi th other observations with this instrument {e.g. 
[Langangen et al. 2007 ). This relative straylight level is then sub- 
tracted from each spectrogram. 



3.4. Normalization 

Normalization of the spectrogram is one of the most important 
steps of the reduction. There are several factors that determine 
the overall shape of the spectra: filter transmission function, grat- 
ing blaze function, vignetting, fringing, etc.. Flat fielding re- 
moves most of the fringing, however our procedure will leave 
traces of fringe components parallel to the spectral lines. Our 
approach to normalization is a 'one-fits-all', in which we cor- 
rect for all these eff'ects in one step. Once again, the FTS atlas 
is used. It is assumed that it perfectly reproduces the mean spec- 
trum at disk-centre. Our procedure is to find the transformation 
that maps the continuum of our mean spectrum into the contin- 
uum of the FTS atlas. This is done by first convolving the at- 
las with a Gaussian instrumental profile. Our mean spectrum is 
then divided by the atlas and a 'continuum map' is obtained by 
smoothing this ratio. We find that this procedure is extremely ef- 
ficient at removing fringing and other artifacts. Due to the filter 
position the wavelength region around 630 nm happened to be 
the most aff'ected by fringing. Figure [2] gives an example of how 
the continuum finding procedure works in this region. 

Because there is no absolute intensity calibration, normaliza- 
tion of the calibration mean spectrum defines the disk-centre in- 
tensity. All spectrograms are divided by the corresponding con- 
tinuum map, which traces the continuum level in the mean spec- 
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Fig. 2. Normalization of the mean spectrum, for a frame in the 
630 nm region aff'ected by fringing. Top: mean spectrum before 
normalization (thin line), continuum map from our algorithm 
(thick line). Bottom: mean spectrum after normalization. 



trum. This ensures that all the distortions in the spectra are cor- 
rected. The next step is to find the continuum level, and for each 
spatial point, as it varies along the granulation pattern. This is 
done by identifying the continuum regions in the spectrum (de- 
fined as points where / > 0.995 in the FTS atlas). The continuum 
level is taken as the mean intensity in such regions. 



3.5. Frame selection and sampling 

In our observations more than 392 000 spectrograms were ob- 
tained. From these, only about 500 were actually used in the 
analysis. The most important criterion for data selection was the 
continuum intensity contrast. The maximum values for the con- 
tinuum contrast for the 615 nm are shown in Table [T] The disk- 
centre value of 9.2% demonstrates the excellent quality of the 
data. It can be compared with other values obtained with slit 
spectrographs, like the ^ 7% at 615 nm of Kiselman (1994 ) with 
the SV ST and the 7% at 630 nnj^ obtained by Da nilovic eraT] 
^OE^ with the SOT/SP of the Hinode space telescope. 

Continuum contrast is a good measure of quality for quiet 
Sun regions, but cannot be used blindly. Some magnetic fea- 
tures are difficult to identify just by visual inspection of the slit- 
jaw images, and can artiffcially increase the continuum contrast. 
Because the present work is concerned with quiet regions, some 
subsets had to be removed from the analysis because they were 
suspect of having magnetic eff'ects (e.g., excessively bright spa- 
tial points). Other sets displayed an excessive amount of fring- 
ing that could not be completely corrected, and therefore were 
discarded. The best subset for each yu value was used. From each 
subset selected for analysis the best 25 images were used, except 
for disk-centre where the best 50 images were used. 

In Figure |3] an example spectrogram is displayed, showing 
the eff'ects of reduction in the image. 

An important aspect of the selected frames is the temporal 
and spatial sampling. With the simulations we have a very good 
time and spatial sampling of the solar granulation. However this 



^ The continuum contrast is inversely proportional to wavelength (ap- 
proximately). Our value of 9.2% at 615 nm corresponds to ^9.0% at 630 
nm. 
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Fig. 3. Example spectrogram for the 777 nm region before (top) 
and after (bottom) the reduction. In the final reduced image dis- 
tortions and inhomogeneities such as the fringing pattern in the 
top right; vignetting; continuum function and dust specs in the 
CCD have been corrected. The vertical axis represents spatial 
position and the horizontal axis wavelength. The three vertical 
lines on the left are the O i 777 nm lines, and the darker line on 
the right is the strong Fe i 778.05 nm line. 



is not always possible with the observations if we want to se- 
lect the best continuum contrast. The observations with higher 
contrast will tend to be clustered in the instants with better see- 
ing, which do not provide the best sampling. There is a trade-off' 
between getting the highest quality images and getting the best 
sampling. In this work we have seeked the higher quality images, 
at the expense of sampling. 

A spatial sampling as good as the simulations will always 
be difl&cult to achieve because the models have a two dimen- 
sional sampling of the solar surface, while the spectrograph's 
slit is only one dimensional. Our slit covers ^ 40 arcsec, while 
the models cover a box of 8 x 8 arcsec^. In the time domain 
the simulations have a better coverage (20 snapshots uniformly 
sampled from 45 minutes of solar time). Our observations at 
disk-centre cover slightly more than 10 minutes, while at other 
viewing angles the time coverage is typically of a few minutes. 
The solar 5 -minute oscillations will inffuence the line properties 
(mainly velocities), so a good temporal sampling is important to 
average the oscillations, which are also present in the simulation. 
For our disk-centre data the coverage is adequate for this pur- 
pose, but for other viewing angles it is not so good. The problem 
of systematics in the line shifts for yu 9^ 1 due to solar oscillations 
has been minimized by measuring the line shifts relative to the 
spatially and temporally averaged spectrum for the given yu, and 
not from the mean disk-centre spectrum. 




615.40 615.45 615.50 615.55 615.60 615.65 615.70 615.75 615.80 615.85 
Wavelength [nm] 



Fig. 4. Eff'ect of the Fourier filter in a spectral region around 615 
nm. The high-frequency noise components are eliminated in the 
filtered spectrum. 



3.6. Fourier filtering 

The signal-to-noise ratio of each individual spatially-resolved 
spectrum is on the order of 25-50, depending on the camera 
used and the exposure level. To minimize this noise a Fourier 
filter was implemented. Because the spectrograms are oversam- 
plecj^in the spatial and wavelength directions a Fourier filter is 
very eff'ective at minimizing the noise with little loss of signal. 
Given that on the spectrograms noise has a much higher spatial 
frequency than the spectral lines, a low-pass filter eliminates the 
noise while preserving the spectral information. Working on the 
Fourier space the high frequencies are filtered out by multiplying 
the two dimensional Fourier transform of the spectrogram by a 
low-pass filter function, and taking the inverse Fourier transform 
to obtain the filtered data. For the filter function we use a two di- 
mensional Butterworth filter with an elliptic shape, defined as: 



(3) 
(4) 



where n is the order of the filter, D is the distance in the two 
dimensional frequency space (fx, fy), and Cy are the cutoff* 
frequencies in x and y. An order of = 20 was used, and the 
chosen cutoff* frequencies were c^^ Cy = O.lpix"^ These values 
ensured that a significant part of the noise was eliminated while 
retaining the spectral information. The smooth dip in the func- 
tion ensures that Fourier artifacts such as the Gibbs phenomenon 
are minimized. Being a two dimensional problem, noise in both 
axis of the spectrogram (wavelength and spatial coordinate) can 
be treated diff'erently by setting diff'erent values for Cx and Cy. 

Even with the conservatively chosen cutoff' frequencies, 
Fourier filtering significantly increases the signal-to-noise ratio. 
Its eff'ects are very helpful for a precise continuum determination 
and for very weak lines. Fig. [4] shows an example of the eff'ects 
of the filtering. 



^ There are at least 5 pixels per resolution element in the spatial axis 
and 3 in the wavelength axis. 
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3. 7. Extracting line quantities 

Five line properties are extracted from the spectra: line-centre 
position, line-centre intensity, full width at half maximum 
(FWHM), a measure of line asymmetry, equivalent width. These 
are computed automatically (i.e., no human input in the fitting 
procedures) using algorithms designed for the eff'ect. 

We define 'line centre' as the velocity diff'erence in kms"^ 
between the line centre of the mean spectrum (spatial and tem- 
poral average the spectra at the given ji value) and the line centre 
for the spectrum of a given spatial point. The line centre wave- 
length and intensity are obtained by fitting a parabola to a few 
points around the line minimum. The velocity shift is then com- 
puted using the mean spectrum and converting wavelength to ve- 
locity units. The line centre intensity is given in relative intensity 
units. 

The line FWHM is computed by fitting the line wings and 
measuring the width at half maximum. Two 5th order polynomi- 
als are fitted to the red and blue wings of the line being analyzed. 
The fitted functions are then evaluated at the half maximum 
of the line, and the wavelength distance between red and blue 
wings is measured. This distance corresponds to the FWHM, 
and is converted to velocity units. The quantity 'line asymme- 
try' is computed in a similar way. The line asymmetry is defined 
as diff'erence (in velocity units) between the line centre and the 
bisector at half maximum. The wavelength position of the line 
bisector at half maximum is obtained as the FWHM: from the 
fitted polynomials to the line wings. It is then subtracted from 
the line centre wavelength and converted to velocity units. Due 
to the fitting procedure, the FWHM and 'line asymmetry' are the 
quantities with the highest uncertainties. 

The equivalent widths were computed by direct integration. 
For consistency both the observations and the (convolved) sim- 
ulations were run through the same analysis procedure. 

The uncertainties for all the line quantities reflect the pho- 
ton noise and continuum placement errors. The efl'ects of the 
photon noise, were estimated using artificial lines. To generate 
an artificial line with similar properties, a Gaussian was fitted 
to the mean spectrum of each line. This artificial line was then 
placed in a continuum region in a few spectrograms. Its quan- 
tities were computed in the usual fashion and the comparison 
between the Gaussian exact quantities with the extracted quanti- 
ties yields the errors due to photon noise. The final errors were 
the quadratic sum of the photon noise and continuum placement 
errors. The continuum placement error was assumed to be 0.25% 
for all lines, for each spatially-resolved spectrum. 



4. Simulations 

4.1. Atmospheric model 

The 3D line formation calculations presented herein have been 
performed with a 3D, time-dependent, radiative-hydrodynamical 
simulation of the solar photosphere (Trampedach et al., in prepa- 
ration; Asplund et al., in preparation). This new 3D mo del differs 
in a number of ways from the 3D solar atmosphere of 'Asplund' 
|et al.| pQQQ) that has been extensive ly used for solar abun dance 
purposes over the past decad e {e.g. AUende Prieto et aL||2004] 
[Asplund et al.||20Q4[ |2005b|). It has been computed with the 



(vertical direction plus 2 inclined with 4 azimuthal rays). The 
code makes use of periodic horizontal and open vertical bound- 
aries. The Eulerian mesh consists of 240^ cells, which corre- 
sponds to a physical dimension of 6x6x4 Mm; the solar simula- 
tion extends about 900 km above the optical depth unity. For the 
subsequent line formation calculations the original simulation 
was interpolated to a finer vertical depth scale while sampling 
the horizontal directions to produce a 3D model of dimensions 
50x50x82; test calculations confirmed that this procedure did 
not introduce any diff'erences in the final results. 

Compared to our previous 3D solar model, an extensive over- 
haul of the input micro-physics and radiative transfer treatment 
has been carried ou t. The continuum and l ine opacities are based 
on the MARCS code ( [Gustafsson et al. 2008 ) supplemented by ad- 
ditional photo-ionization cross-sections from the Opacity and 
Iron Projects ( |Cunto et aL||1993[ [Hummer et~ar| 11993b . The 
adopted equation-of- s tate is that of jMihalas et al. jl988 ). The 
[Asplund et al. ( 2005a I solar chemical composition has been used 
both for the opacities and equation-of- state. The wavelength- 
dependence of the opacities have be en accounted for using the 
concept of opacity binning (jNordlund 1982). The division in 
opacity and wavelength for the 12 opacity bins was based on 
the full monochromatic (105 wavelength points) solution for the 
mean 3D stratification. Continuum scattering was included as a 
pure absorption (i.e. S ^ = B;i(T), which is fully justified for the 
solar photosphere (Hayek et al., in preparation). 

The time- sequence of the 3D solar simulation utilized for the 
line formation calculations corresponds to 1.2 hr of solar time. 
The resulting eff'ective temperature of the 90 snapshots is Tgff = 
5778 + 14 K. 

4.2. Line formation and departures from LTE 

Using the 3D model described above we computed the synthetic 
line profiles using our LTE line formation code. From the full 
simulation a series of 20 snapshots covering ^40 min of solar 
time has been used. The adopted snapshots were chosen so that 
the temporally averaged eff'ective temperature of the atmosphere 
(not a free parameter in the simulation) was Tqq = 5777 + 14 K, 
very close to the observed nominal value of 5777 ± 3 K ( [Willson 



STAGGER MHD c ode ( [Nordlund & Galsgaard|1995[[Gudiksen &[ 
|Nord lund" 2005 ), although no magnetic field has been imposed 
for the present simulation. The code solves the hydrodynamical 
equations of mass, momentum and energy coupled with a solu- 
tion of the 3D radiative transfer equation using nine directions 



& H udson, 1988| ). While the full simulation has a numerical res- 
olution of 150 X 150x 82, for the purposes of spectral line forma- 
tion we have used an interpolated lower resolution of 50x50x82. 
As outlined in [Asplund et al.[ ( [2000) this procedure is meant to 
save computing time and does not introduce any diff'erences in 
the line profile shapes. Assuming LTE, the line formation code 
used a Planckian source function (Sy = By), and the Saha and 
Boltzmann equations to determine respectively the ionization 
fraction and atomic level populations for each species. For realis- 
tic equation-of-state and continuum opacities the Uppsala opac- 
ity package was employed. The radiative transfer equation has 
been solved for several directions: for each of the five ji posi- 
tions of our observations the line profiles were computed over 
for four (f angles. For each line the profiles were computed us- 
ing three abundance values, allowing for interpolation of other 
abundance values. No additional line broadening through micro 
or macroturbulence was applied. 

The synthetic line profiles were computed assuming the LTE 
approximation. For oxygen lines it was possible to perform 3D 
NLTE line formation calculations. This was accomplished us- 
ing the multi3d cod e ( Botnen[l997t Botnen & Carls son][l 999 [ 
[Asplund et al.|2003] ), which solves the rate equations iteratively 
using a local lambda operator and assuming statistical equilib- 
rium. For the purposes of calculating the level populations, the 
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radiative transfer in multiSd was solved for 40 outgoing inclined 
rays using a short characteristic technique. Once the level pop- 
ulations had converged and the iteration is finished, we made 
use of the resulting source function and opacities to compute 
the radiative transfer using a long characteristic Feautrier solver. 
Similarly to the LTE code, we compute the NLTE line profiles 
for the [I values of our observations, with 4 (f angles each. The 
Uppsala opacity package was used to compute the continuum 
opacities and equation-of- state for multiSd. For increased nu- 
merical accuracy, the 3D model fed into multiSd had a resolu- 
tion of 100 points in the vertical scale, interpolated from the 82 
points of the original 3D simulation. 

Due to the computational expense, full 3D NLTE line for- 
mation was carried only for 10 snapshots of the simulation. The 
snapshots were selected to provide a good statistical coverage, 
both in time (evolution of granulation) and eff'ective tempera- 
ture. The final NLTE line profiles were obtained by multiply- 
ing the LTE line profiles (computed with our LTE code) by 
the NLTE/LTE ratio (computed with multi3d). The wavelength- 
dependent NLTE/LTE ratio was obtained for each snapshot, 
viewing angle and spatial point. Multiplication of the ratio by 
the line profiles is made individually for all spatially-resolved 
spectra and for each viewing angle and snapshot. 

Because one has NLTE profiles directly from multi3d, it may 
seem counterproductive to obtain the final values by multiplying 
the LTE profiles by the NLTE/LTE ratio. However, this approach 
was followed for two reasons. First, to allow interpolation of the 
line profiles in abundance. Line profiles from multi3d are com- 
puted only for one abundance value, whereas for our LTE code 
they are computed for three values of log(6o • gf). Interpolating 
between these three line profiles allowed us to calculate the line 
profile for an arbitrary oxygen abundance. The second reason 
this approach was taken was to ensure consistency between dif- 
ferent lines by using the same line formation code. In any case, 
for the oxygen lines studied we find very little diff'erence in the 
LTE line profiles between multi3d and our LTE code. 

4.3. Atomic data 

The two main sources of a tomic data for our spec tral calculations 
were the NIST database (jRalchenko et al.|2009l) and t he VALD 
database ( [Piskunov et al.||1995[ |Kupka et al.||2000| ). A sum- 
mary of the lines properties and the sources used can be found 
in Table [2] CoUisional (van der Waals) bro adening was com- 
puted using the quantum mechanica l theory of|Anstee & O'Mara 



T995l ); [B^lem & Q'Maral ( [T997l );'Barklem et al.|( |1998| ), elim- 
inating the need for conventional [Unsold (195 51 enhancement 
factors. The collisional broadening coefficients were computed 
for the oxygen lines and, when available, to other lines used in 
our analysis. 

Aside from the lines listed in Table [2] we have included 
nearby blends in some line profiles. For the [Oi] 630.03 nm 
we included the Nil blend and for Oi 615.81 nm we have in- 
cluded some weak molecular and atomic blends (f or more de- 
tails about the blends in these lines, see [Paper II| ). In the line 
profiles of Si I 615.57 nm and Cai 616.60 nm we have included 
both lines and also a nearby strong Sii line at 615.51 nm. In 
the case of Cai 616.60 nm we have also included the very weak 
Oi triplet at 615.60 nm. Data for all these lines were extracted 
from the VALD database. For the Nil 630.03 nm line we use 
\oggf = -2.11 from Johansson et al.|( [20Q3] ). 

For the NLTE calculations we employ a 23 level O atom 
(22 O I levels and the ground state of O ii), including in total 43 
bound-bound and 22 bound- free transitions. This atom is based 



Table 2. Lines studied in the present work and their parame- 
ters. Equivalent widths at disk-centre measured from our obser- 
vations. Data sour ces : a) VALD database; b) N IST database; c) 
Nave etal.| ( |T994l ); d) |Storey & Zeippen| ( |2Q5Ql ). 



Atomic Species 


A [nm] 




^low [eV] 


W,=i [pm] 


Fei 


615.1618'^ 


-3.299'^ 


2.18'^ 


4.95 


Si I 


615.5693^ 


-2.252^ 


5.62" 


0.61 


Cai 


1 c /^rvooa 

oij.bUzj 


-ZAy / 


Z.JZ 




Fei 


615.7728'^ 


-1.260" 


4.08^ 


6.45 


Oi 


615.8186^ 


-1.841^ 


10.74^ 


0.63 


Fei 


615.9378'^ 


-1.970" 


4.61" 


1.21 


Fei 


629.0965^ 


-0.774" 


4.73" 


7.06 


Fei 


629.7792^ 


-2.740^ 


2.22^ 


7.54 


[Oi] 


630.0304^ 


-9.717^ 


0.00^ 


0.45 


Sen 


630.0698^ 


-1.887" 


1.51" 


0.59 


Oi 


777.1944'^ 


0.369'^ 


9.15'^ 


8.48 


Oi 


777.4166^ 


0.223^ 


9.15^ 


7.44 


Oi 


777.5390^ 


0.002^ 


9.15^ 


5.90 


Fei 


778.0557^ 


0.029" 


4.47" 


14.81 



on the atom of Kiselman (1993]), but has undergone updates 
( Fabbian et al.|200 9) as new data for electron-impact excitation 
has recently become available (Barklem^2007j . 

4.4. Image degradation and instrumental profiles 

To compare the simulations with Earth-based observations one 
needs to degrade them in a similar manner as the Earth's 
atmosphere and the telescope degrade the images from the 
Sun. Additionally, the spectra will also have to be con- 
volved to account for the spectrograph's instrumental pro- 
file. We employ a three step convolution process that mimics 
the path of light from the solar surface to the spectrograph: 
atmosphere telescope spectrograph. The first step is the 
point-spread function (PSF) caused by the Earth's atmosphere. 
Following Nordlund| ( [T984| ) and jCoUados & Vazquez] fl 987) we 
describe the atmospheric PSF as a sum of two Lorentzians: 



/?atm(^ a, b) = 



(5) 



where r is angular distance and we define a > b so that the first 
Lorentzian determines the overall contrast ('extended wings' of 
the PSF) and the second Lorentzian determines the size of the 
smallest detail observable ('narrow core' of the PSF). Tests using 
diff'erent functions for the PSF show that the eff'ect on the results 
is negligible. 

Following the path of light, the next step is to mimic the finite 
resolution of the telescope. We employ an Airy disk function for 
a circular aperture telescope as the 'diff'raction' PSF: 



D 



(6) 
(7) 



where D is the telescope's diameter, A the wavelength of the ob- 
servations, r the angular distance and Ji a Bessel function of the 
first kind. 

The combined PSF is then the convolution of p^tm with /^diff : 



p(r, a, b) = [/7atm * Pm] (r, a, b). 



(8) 
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Fig. 5. Convolution of the simulations with the telescope's 
and the Earth's atmospheric PSF at the solar disk-centre: a) 
Continuum intensity at 615 nm for the original simulation; b) 
Continuum intensity at 615 nm for the convolved simulation; c) 
The middle row of the 2D PSF used for the convolution (con- 
verted from simulation pixels to arcsec and in log scale); d) nor- 
malized continuum intensity distribution at 615 nm for the sim- 
ulations and observations. The convolution parameters a and b 
(see text) are chosen so that the continuum intensity distribution 
and continuum contrast of the simulations match the observa- 
tions. 



The PSF is normalized, so that the convolution with the im- 
ages conserves the intensity. The a and b parameters depend on 
the local conditions and are empirically obtained. Two properties 
of the observations are used to constrain them: the continuum in- 
tensity contrast and the continuum intensity distribution. These 
two properties are not completely independent of each other: 
a change of the distribution will change the contrast and vice- 
versa. Nevertheless it is possible, for a fixed contrast, to change 
the shape of the distribution using diff'erent pairs of (a,b), and 
that is the approach taken. For the disk-centre data we adjust 
a and b in order to get the best agreement with the observed 
continuum intensity distribution and continuum contrast. Fig. |5] 
shows an example of how 777 nm data is convolved so that it 
agrees with the observations. By fitting the observed intensity 
distributions, the adopted PSF formulation also corrects for the 
non-ideal components of the telescope and instrument PSFs. 

There is some degree of degeneracy in the intensity distri- 
bution for diff'erent values of a and b. The adopted PSF is not 
necessarily the most physically correct. It is rather an empirical 
transformation that best reproduces the observations, given the 
coarse spatial resolution of the simulations and the time-average 
of diff'erent seeing conditions that the observed intensity distri- 
bution represents. Several tests were made to study the eff'ect of 
diff'erent PSFs in the results. In Fig. [6] we plot some results for 
two sets of parameters giving diff'erent PSFs: one with stronger 
wings (left column), and the other with weaker wings (right col- 
umn). The diff'erence in the intensity distributions given by these 






2 - 


ms 


1 - 
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-2 - 



0.8 0.9 1.0 1.1 1.2 0.8 0.9 1.0 1.1 1.2 
Continuum intensity Continuum intensity 

Fig. 6. Eff'ects of diff'erent PSF shapes at 615 nm, ji = 1, us- 
ing two sets of a and b parameters (left and right columns). Top 
panels: middle row of the 2D PSF. Middle panels: continuum in- 
tensity distribution for the convolved simulations (thin line) and 
the observations (thick line). Bottom panels: eff'ects of the diff'er- 
ent convolutions in the line shifts of the strong Fei 615.77 nm 
line. Contours represent the point density (see Sect. 5.1 ). Solid 
colour contours are for the convolution with the PSF adopted for 
our analysis (see Fig. [5]), line contours for the convolution with 
the PSF of the corresponding parameter set. 



two PSFs is very small, as is the eff'ect on the results (in this case, 
plotted for the line shifts, but the results are similar for other pa- 
rameters). We conclude that the adopted PSF describes the ob- 
servations reasonably, and the error introduced in the results by 
a diff'erent shape of PSF is small. 



Stein & Nordlund ( 1998 ) hint at a good agreement between 
the 3D model intensity distibutions and the observations. Recent 
work has indeed shown that the intensity distribution of these 
and other simulations are consistent with observations once care- 



ful PSF modelling is performed (Wedemeyer-Bohm & Rouppe 
van der Voort|20Q9[|Danilovic et al.|2Q08| ). |Wedemeyer-Bohm & 
Rouppe van der Voort| ( |2Q09 l) also show that the observed centre- 
to-limb variation of the continuum contrast is well described by 
a 3D model of the Stein & Nordlund family.. Nevertheless, and 
for a fair comparison, we chose not to assume too much about 
the predicted intensity distributions at yu 9^ 1. Unlike at disk- 
centre, we do not force the model intensity distribution to match 
the observations. The parameter a is related to the PSF wings. 
Its sources are several: stray light in the telescope, atmospheric 
scattering, etc. It is reasonable to assume that it will not change 
significantly between observing sets. The parameter b, on the 
other hand, defines the seeing-induced resolution, which is ex- 
pected to vary between sets. Hence for the yu 9^ 1 sets we use 
the same a value as obtained for disk-centre, and adjust only b 
to obtain the observed continuum contrast. 
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The simulation used has a horizontal resolution of 50 x 50. 
For disk-centre this corresponds to a spatial scale of 6 x 6 Mm. 
Because the pixel size of the box is kept constant, when one com- 
putes line profiles at yu 9^ 1 the spatial resolution will increase 
with decreasing ji. Diff'erent ip angles at off'-centre positions will 
also cause diff'erent spatial resolutions along diff'erent axes of 
the simulation box. The two-dimensional PSF is adjusted so that 
the angular coordinate r matches the spatial coordinate in pix- 
els, meaning it has an ellipsoidal profile to account for yu, rotated 
according to the (p angle. 

After convolving the synthetic spectrograms with the PSF, 
each individual spectrum is also convolved with the instrumental 
profile of the spectrograph. A Gaussian equivalent to a resolving 
power of 200000 was used. This value of the resolution was 
obtained by investigating the data and comparing with the FTS 
intensity atlas. 



5. Results and discussion 

5.1. Line shifts 

The line shifts are a diagnostic of the velocity fields on the solar 
surface. Observations at disk-centre probe the vertical motions, 
while observations closer to the limb can also provide some in- 
formation about horizontal motions. It should be noted, however, 
that because line formation takes place at a range of depths it is 
not possible to associate the observed line shifts with a specific 
geometric height. The lower contrast of the off'-centre data (espe- 
cially at the limb) means that a precise comparison of the line- 
of-sight velocities between the diff'erent positions on the solar 
disk is very difficult. 

Some results are plotted in Fig. |7] The diff'erent lines have 
diff'erent formation ranges, and the line shifts give us informa- 
tion about the velocity fields where these lines were formed. 
For the line shifts as well as most of the other line information, 
the disk-centre data provide the most detail in the structure as 
a function of continuum intensity. One can see very clearly the 
trend with the continuum intensity: the dark (intergranular) re- 
gions representing downflows and the bright (granules) regions 
representing upflows. The agreement with the 3D model is very 
good, confirming its realistic velocity fields. 

The results for other lines are very similar to Fig. [7] with 
generally an excellent agreement between the 3D model and the 
observations. For some very weak lines with close stronger lines 
(e.g. 1 615.81 nm) our algorithm for determining the line centre 
can fail more often. This is due to the low S/N ratio and difficulty 
in finding a very weak line, and results in an increased scatter in 
the line shifts. 



5.2. Line strengths 

The distribution of the line equivalent widths over the solar gran- 
ulation is one of the most interesting diagnostics. It provides 
a wealth of information on how lines are formed over a depth 
range. From an observational point of view it can be used to 
'finge rprint' particular types of species (K iselman & Asplund 
2003| ). But more interestingly, it can be used to probe the struc- 
ture of 3D models and the microphysics of line formation. The 
early work of Kiselman ( |1998| ) is an example of how similar ob- 
servations can be used to probe for d epartures from LTE in the 
671 nm Lii line. Indeed in Sect. 5.5 we detail a similar analysis 
of the NLTE eff'ects of the O i lines using the spatially resolved 
variation of equivalent widths. 



In Fig. [8] we show the distribution of equivalent widths for 
all the lines considered in the present study. The contours reflect 
the density of points in the plot. Statistically observations and 3D 
model are comparable, the contours comprising ^^50 000 spatial 
points, exception made for the O i 777 nm lines, where the con- 
tours from the NLTE profiles comprise 25 000 points (only 10 
snapshots were used for NLTE calculations). 

With a relatively low signal-to-noise ratio for each spatial 
point, it is not surprising that the observations of weak lines 
show a bigger scatter than the 3D model. Photon noise influ- 
ences the equivalent widths directly and also the determination 
of the continuum level for each spatial point, efl'ects that become 
more pronounced in weaker lines. 

Overall the agreement between the 3D model and the ob- 
servations is very good. The Sii 615.57 nm line shows a slight 
mismatch, with the model predicting a weaker line at the bright 
end. This may be due to blend contamination. In the line profile 
we have included blends with a nearby strong Si i line and also 
the Ca I 615.60 nm (because their wings overlap), but the prob- 
lem persists. In this region the VALD database lists only a very 
weak Fe i which is unlikely to alleviate the problem either. We 
note that the Ca 1 615.60 nm would display a similar behaviour if 
we had not included the blend with the O i 615.6 nm line (£'iow = 
10.74 eV). For the Si i it is possible that a similar high-excitation 
blend might be at play. 

For the other lines it is interesting to note the eff'ect of the 
excitation energy: the 1 777 and 615.82 nm having a high exci- 
tation energy and a very pronounced equivalent width variation, 
while other lines have a more 'C shaped pattern, becoming in- 
creasingly flat as the excitation energy decreases (e.g. [O i]). The 
variation of the equivalent width distribution with excitation en- 
ergy is discussed in Sect. |5.6[ For the strong Fe i lines the 3D 
model seems to predict a slightly steeper variation of equivalent 
width with continuum intensity at the bright end. This is visible 
for the Fei 615.77 nm, 629.09 nm, 629.77 nm and, to a lesser 
extent, 778.05 nm. The origin of this efl'ect is unclear. 

The line profiles were adjusted in abundance so that they 
match the observed Wmspec- This adjustment essentially moves 
the distribution vertically in Fig. [8] Again, it is very reassuring 
to see how well the 3D model and LTE line formation (except 
O I 777 nm) describe the observed distributions. 

5.3. Line asymmetries 

The line-asymmetry measure is defined here as the wavelength 
diff'erence (in velocity units) between the line core and the bi- 
sector at half maximum. It is a two-point sampling of the line 
bisector and, along with the line shift, another probe of the ve- 
locity fields. For a given point in the solar granulation, the line 
asymmetry reflects the velocity distribution over its formation 
depth and optical path. 

In Fig. [9] we show the variation of line asymmetry at disk- 
centre for four medium to strong lines. For weak lines the asym- 
metry distribution is essentially flat, meaning that the velocity 
difl'erences along the formation region are small. The Oi 777 
nm lines provide a great probe for the line asymmetries, be- 
cause they have essentially the same line parameters except the 
log gf value. In Fig. [9] we can see how the asymmetry distri- 
bution changes from the the weakest line (777.53 nm) to the 
strongest (777.19 nm). For the weakest line there is almost no 
variation with continuum intensity, and an overall positive asym- 
metry. This means that the 'bulk' of the line profile is blueshifted 
in regards to the line core, a result of the vertical velocities in its 
formation region being mostly positive (outgoing material). As 
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Fig. 7. Line shifts for selected lines and several jd values. Data shown as contour levels of the two dimensional histogram of the 
density of points. Observations are depicted as solid colour contours, while the 3D model as line contours. The continuum intensity 
is normalized to the mean intensity for each ji value. The lower continuum contrast towards the limb is evident, with less structure 
information at yu = 0.2. The energy of the lower level is shown, along with the disk-centre equivalent width. The chosen lines cover a 
representative range of lower level energies and equivalent widths. The results for O i 615.81 nm are uncertain due to the line being 
particularly weak, especially at intergranular regions and at the limb. Results for the other two O i 777 nm triplet lines are similar 
those shown for O i 7774.1 nm. 



the lines get stronger their formation range will extend outwards. 
The vertical velocities in higher layers of the photosphere tend 
to get progressively smaller (with an increase of the horizon- 
tal motions as the gas outflows cool down). Formed at higher 



depths, the stronger lines will cover a larger velocity gradient in 
the photosphere, and a trend of asymmetry with continuum in- 
tensity becomes noticeable. The dark intergranular regions will 
have a negative asymmetry (the line profile is redshifted in re- 




Fig. 8. Distribution of equivalent widths over the solar granulation at disk-centre. For each spatially-resolved spectrum the equivalent 
width has been computed and the resulting histogram of points is represented by the contours, as a function of the normalized local 
continuum intensity. The line profiles of Sii 615.57 nm, Cai 615.60 nm, Oi 615.82 nm and [Oi] 630.03 nm include nearby blends 
(see Sect. |4.3| ). All line profiles were computed assuming LTE, except the O i 777 nm lines, showing the NLTE results for 5'h = 1. 
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Fig. 9. Line asymmetries at disk-centre for four of the strongest 
lines studied: the Oi 777 nm triplet and Fei 778.05 nm. 
Comparing observations (solid colour contours) with 3D model 
(line contours). For Oi 777 nm showing results from LTE pro- 
files due to higher number of snapshots computed (the asymme- 
try results are not very sensitive to departures from LTE). 
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Fig. 10. Line asymmetries vs. line shifts for the 3D model, LTE, 
disk-centre and convolved to match the observations. For clarity 
plotted with points instead of contours. Dark points represent 
/c > (/c), light points k < (h)- 



lation to the core) and the brighter regions will have a positive 
asymmetry (the line profile is blueshifted in relation to the core). 
These eff'ects are a result of the downflows and upflows of gas 
at these regions. They become even more apparent in the much 
stronger Fei 778.05 nm line, covering a significant velocity gra- 
dient and displaying the most extreme values of asymmetry from 
all lines. 

For these four lines as well as most of the other lines (not 
shown) the 3D model predictions regarding line asymmetry are 
very good. This is not surprising, given the excellent results re- 
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Fig. 11. Line full widths at half maximum (FWHM) at disk- 
centre for four lines of different strengths and excitation ener- 
gies. Comparing observations (solid colour contours) with 3D 
model (line contours). 



ported earlier about the line bisectors ( [Asplund et al.|2000| ). Our 
results for line asymmetries are consistent with the findings of 
|Kiselman |([T994). 

In Fig. [To] a diff'erent perspective is shown. Here the asym- 
metries of the weakest of the 1 777 nm lines and the strong Fe i 
line are plotted against the line shifts (measured from the line 
core in relation to the mean spectrum). The points are colour- 
coded between dark and bright continuum intensity regions. For 
the weaker line the average line shift for the bright granules is 
negative, and positive for intergranular regions. For the strong 
line granules and intergranular lane regions occupy roughly the 
same space in line shift, and the overall dispersion in line shifts 
is smaller. One explanation for this could be that the core of 
this line is formed above the convective overshoot region. But 
while the core is formed higher, the asymmetry of the line pro- 
files reflects the motions below, and it becomes very clear that 
the negative asymmetry is associated with dark regions and pos- 
itive asymmetry with bright regions. This distinction is hardly 
visible in the weaker line, formed deeper in the photosphere. 

5.4. Line FWHM 

The distribution of the FWHM with continuum intensity for 
disk-centre is shown in Fig. [TT] for four lines. For most of the 
lines the FWHM decreases with increasing continuum intensity, 
with a larger scatter at low continuum intensity. This decreas- 
ing behaviour of FWHM with continuum intensity has been ob- 
Hanslmeier eFaLlp^93l |Kiselman|[T994l 



served by others (e.g. 



Puschmann et al.|2QQ3| ) and is caused by increased line broaden- 
ing due to turbulence in the intergranular regions. For example, 
tests with the [O i] 630 nm line show that if the velocity fields 
are switched off' when computing the theoretical profiles, then 
the predicted FWHM distribution is fiat. Very strong lines like 
Fei 778.05 nm are formed higher in the atmosphere and above 
the convective overshoot region, leading to a smaller turbulent 
broadening and a ffatter FWHM distribution, as seen in Fig.[TT] 



The O I 777 nm lines are a special case in the FWHM distri- 
bution. Unlike all the other lines considered here, their net varia- 
tion of FWHM is an increase with the continuum intensity. It can 
be seen in the bottom panel of Fig. 12 for Oi 777.41 nm. This 
eff'ect is related to thermal broadening and results from a pecu- 
liar combination of strong lines with a high excitation potential 
and the relatively light atom of oxygen. With £^iow = 9.15 eV the 
Oi 777 nm lines are formed very deep where the temperature 
is higher. A high temperature together with the relatively small 
atomic weight of oxygen causes the large thermal broadening 
seen in these lines, which dominates over turbulent broadening. 
With a steeper temperature gradient in the granules the thermal 
broadening is larger in the bright regions, driving the FWHM in- 
crease with continuum intensity. Tests with fake lines show that 
the negative slope of the FWHM with continuum intensity can 
be recovered if the atomic weight is assumed to be much higher 
than that of oxygen (thereby minimizing thermal broadening), 
confirming that the FWHM trend is dominated by thermal broad- 
ening. 

Another candidate for this kind of behaviour is the Oi 
615.81 nm (E\o^ = 10.74 eV). However, because this line is 
much weaker, in absolute terms its thermal broadening will be 
much smaller, and turbulent broadening dominates. The model 
predictions (not shown) confirm a decrease of FWHM with 
the continuum intensity. The inclusion of the molecular blends 
around this line (see discussion in Paper II ) will reinforce even 
more the negative variation of FWHM with continuum intensity. 
In the observations the uncertainty in its FWHM is very large 
because the line is very weak, preventing a confrontation with 
the model predictions. 

Overall the 3D model reproduces the observed distributions 
of FWHMs very well, for both strong and weak lines and includ- 
ing the special case of the O i 777 nm lines. 

5.5. Departures from LTE in oxygen lines 

Of the Oi lines the 777 nm triplet lines are know to suf- 
fer significant departures from LTE. This was first shown by 
Altrock (1968 ) and later explained by lE nks son & Toft| ( p^79 ). 
Subsequent studies (e.g. |Kiselman] ri991 ) found also NLTE ef- 
fects for the Oi 615.81 nm line, albeit on a much smaller 
scal e - usuall y < 0.03 dex ( [Asplund et"aLl|2004t [Caff'au et aL 
|2008|Paper II| . 

A proper NLTE treatment of these lines poses several chal- 
lenges. As noted in previous studies, one of the main problems is 
to account for the eff'ect of collisions with neutral hydrogen. No 
laboratory data or quantum mechanical calculations are avail- 
able to date. The best estimate to account for the efficiency of 
hydrogen collisio ns has been to use the recipe of Steenbock &| 



Holweger](|1984|), a generalization of the classical Drawin for- 
mula ( |Drawin) 1968 ). This recipe is often scaled by an empirical 
factor 5 H- In the work of|Allen de Prieto et ar| ( |2004) and Pa per II| 



the centre-to-limb variations of the O i 777 nm lines have been 
used to derive the best-fftting Su value. 

We propose to use the line variations across the solar gran- 
ulation at disk-centre to test the LTE and NLTE line formation 
and derive the best fftting empirical Su value. The inherent scat- 
ter of the spatially-resolved line quantities makes this method 
less sensitive or constraining than the centre-to-limb variations. 
Nevertheless, it provides an independent test of the possible Su 
values. 

We have used the same LTE and NLTE calculations as in 



Paper II[ detailed in Sect. |4.2| Due to computational constraints 
the full 3D NLTE line formation was only performed for 10 
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snapshots of the 3D simulation. The NLTE effects in Oi 615.81 
nm are so small that the distributions of line quantities in NLTE 
are essentially indistinguishable from LTE. Hence we limit our 
NLTE analysis to the 777 nm triplet only. The line profiles for 
diff'erent values of 5'h have all been interpolated in abundance 
so that the equivalent width of the spatially and temporally aver- 
aged spectrum matches the observed. 



In Fig. 12 we show how the line equivalent widths and the 
FWHM vary across the solar granulation at disk-centre, for the 
observations and the 3D model with diff'erent Su values for hy- 
drogen collisions. The results shown are for the Oi 777.41 nm 
line, but they are very similar for the other two lines. 

The equivalent width plot is perhaps the most interesting. 
The model results can be divided in two regions: granules and 
intergranular regions, with a smaller point density in the middle, 
where /c = (/c). The hydrogen collisions have diff'erent eff'ects 
on each region. An increase in 5'h makes the equivalent width 
slope less steep in the intergranular regions, but steeper in the 
granular regions. In the LTE case it can be seen that in the inter- 
granular regions the equivalent width slope is steeper than in the 
observations. In such regions the agreement with the observa- 
tions improves for the NLTE models. On the other hand, as 5'h 
decreases the equivalent width slope in the bright regions gets 
steeper, departing from the observed values. No single value of 
Su can be in perfect agreement with the observations in both 
high and low intensity regions. 

Another aspect worth noting is the scatter, especially in the 
dark regions. For Su = 0.01 one can see that the 3D model 
displays a significant amount of scatter at low continuum inten- 
sity. This does not seem to be consistent with the observations. 
As 5'h increases the amount of scatter decreases. This increase 
of scatter in NLTE is the opposite behaviour of the Li i 67 1 nm 
line: for this line there is a much bigger scatter for the LTE pro- 
files ( Kiselman||1998[ [Asplund||2005| ). Because the temperature 
dependence of the Oi 777 nm lines is the opposite of Lii 671 
nm, the opposite scatter results are perhaps not surprising. The 
increased scatter for low 5'h is likely a reflection of the high tem- 
perature sensitivity of the lines. The NLTE efl'ects make the lines 
stronger, and the lower 5'h, the stronger the lines. This strenght- 
ening will increase the temperature sensitivity of the lines, which 
may explain the larger scatter for lower 5'h. If one makes the 
reasonable assumption that the 3D model realistically describes 
dark and bright granulation regions, then the scatter in equivalent 
width is enough to rule out the 5'h ^ 0.3 profiles, when compar- 
ing with the observations. Overall, accounting for the scatter and 
equivalent width slopes in dark and bright regions, 5'h = 1 seems 
to be the value that best reproduces the observations. However, it 
should be noted that for 1 < 5'h < 3 there is not much diff'erence 
in the results. 



In the lower panel of Fig. 12 we show a similar figure but for 
the FWHM. The area occupied by the observations 's contours is 
larger probably due to a better statistical sample. Disregarding 
the line strength, the lower the 5'h, the narrower and deeper the 
NLTE line profile shapes will be when compared with LTE. Thus 
a decreasing Su will drive down the FWHM. This is visible in 
Fig. 12 The LTE FWHM's provide the best match for the obser- 
vations. The NLTE profiles, even for 5'h = 1 are narrower. This 
effect is also seen in the disk-centre profile fits of [Paper n| and 
its cause is unclear. As in the equivalent width figure, lower Su 
show a larger scatter in FWHM. However in this case the scatter 
cannot be used to exclude a particular Sn because the FWHM is 
always lower than the observations. 
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Fig. 12. Top: equivalent widths of the Oi 111 A\ nm line as a 
function of continuum intensity. For clarity the values are shown 
as points, instead of contours. Bottom: FWHM for the same line 
as a function of continuum intensity. Comparing the disk-centre 
observations (dark points in equivalent width; solid color con- 
tours in FWHM) with the 3D model (light points in equivalent 
width, line contours in FWHM), for LTE and NLTE, for different 
values of 5'h. The oxygen abundance has been adjusted in each 
case to fit the observed mean disk-centre equivalent width. 



5.6. Summary quantities of line variations 
5.6.1. Definitions 

In this section we look at the variation of the line quantities at 
high spatial resolution and at several yu values. It is not easy to 
quantify with a few numbers the quantity distributions like in 
Figs.[7]and[8] but it must be done for achieving insight. There are 
many possible correlations between the differ ent line quantities 
can be established. We follow the approach of Kiselman ( |1994 ) 
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Fig. 13. Variation of the a parameter (describing variation of line 
core intensity, see text) with the mean equivalent width at each yu, 
for all lines. The elemental abundances are adjusted so the model 
equivalent widths at yu = 1 match the observed. For yu 9^ 1 the 
disk-centre abundances are used, meaning that the model equiv- 
alent widths depart from the observed in some cases. 



and define the two quantities: 



d/i 



line core 

d/c ' 

(W) ' d/e ' 



(9) 
(10) 



where /une core is the absolute intensity of the line core, W the 
equivalent width and Ic the continuum intensity. Both param- 
eters are obtained by least- squares polynomial fits to the data. 
In the case of a a line is fitted, while for y a parabola is fitted. 
Exceptions to this are the O i 777 nm lines at /i = 1 and yu = 0.8. 
Because their equivalent width variations at these positions are 
not easily described by a parabola, 5th order polynomials are 
fitted instead. 

The a parameter is a measure of the variation of the absolute 
line core intensity with the local continuum intensity. If a = I 
then the line cores follow the continuum variations (expected 
for weak lines). Lower values of a are expected as the lines get 
stronger. This parameter is plotted vs. the equivalent width of 
each line (measured from the spatial mean spectrum at each jj). 

The y parameter describes the variation of equivalent width 
with the local continuum intensity. It is divided by the mean 
equivalent width to put all lines in the same scale. To better de- 
scribe the equivalent width distribution, the derivative dW/d/c 
is evaluated at two points: in the dark (intergranular) regions 
and in the bright (granular) regions. These points are defined as 
the mean continuum intensity of all points where / < (Ic) and 
/ > (/c), respectively. Hence two values of y are defined: ydark 
and rbright- 
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Fig. 14. Variation of the y parameter (describing the variation of 
equivalent width, see text) with the excitation potential for all 
lines, at the diff'erent yu. The parameter is evaluated for dark (top 
panel) and bright regions (bottom panel). 



5.6.2. Results 



It 



The results for a are plotted in Fig. 13 and for y in Fig. 14 
should be noted that these quantities are subject to several un 
certainties (particularly y). The most important source of uncer- 
tainty is the scatter in the quantities and the approximation of 
using a line or parabola to fit these quantities. For t he 3D models 
the particulars of the image degradation (see Sect. |4.4|) also in- 



fluence these quantities, albeit to a lesser extent. The elemental 
abundances have been adjusted so that the model W^^i match 
the observations. The = I abundances are used for all the 
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other fi positions. For the Oi 777 nm lines the NLTE profiles 
with 5'h = 1 were used. 

The a results at yu = 1 are consistent with the findings of 
Kiselman ( |1994j ). The lack of lines with W in the 1-5 pm re- 



gion prevents a clear picture in the variation of a with equivalent 
width, but from the weaker lines a downward linear trend is vis- 
ible. This goes on to a plateau for the stronger lines, where a 
is close to zero or negative. A negative correlation of the line 
core intensity with /c for strong lines has been documented in 
the literature (g.g. [Balthasar et al.|1 990^,'Ha nslmeier et al.|1990t 
[Kucera et al.|1995||Puschmann et al. 2003 ). It is believed to be 
associated with convective overshooting: in higher layers the gas 
is cooler above granules than above the integranular lanes. This 
translates into an inverse granulation pattern, and causes the neg- 
ative correlation between core intensity and /c. The model results 
indicate a slightly lower a for some of the stronger lines, which 
may indicate an higher convective overshooting than observed. 
For fi ^ 1, Of of the strong lines increases. Tests with the 3D 
model at yu = 1 and varying degrees of atmospheric degrada- 
tion show that this eff'ect (and also a less steep dependence of 
a with W) is mostly due to a reduction in continuum contrast. 
Overall the model agrees reasonably well with the observations. 
The stronger lines are more problematic, especially at the limb 
where the predicted equivalent widths diff'er from the observa- 
tions. 

The parameter y measures how the line strength varies with 
the continuum intensity. At = 1 we find both ydark and ybright 
to have a monotonic increase with the excitation potential of the 
lines. This is in agreement with | Kiselman| ( p^94| ), although he 
finds a more clear relation by employing a fil ter to remove the 
eff'ects of solar oscillations. As discussed by |KiselmSi (1994]) 
this trend is an indication of the higher temperature sensitivity 
of lines with higher excitation potentials. A positive y at = 1 
is in agreement with the positive correlation between W and Ic 
found in other studies ( [Hanslmeier et al.|1990t|Puschmann et al.| 
|2Q03j . 

The variation of y between dark and bright regions is small 
but measurable. For most lines ydark > Tbright, indicating a steeper 
variation of W with 1^. The decrease in intensity contrast towards 
the limb (e.g., Fig. [7]) makes the determination of y more prob- 
lematic - the uncertainty from fitting polynomials increases. Yet 
the agreement with the model is good for most lines and yu po- 
sitions. It is interesting to note the agreement between model 
and observations is better for ydark than ybright- The model pre- 
dicts a steeper variation of W with 1^ when 1^ < </c). For the O i 
777 nm lines this is visible in Fig. [121 where for LTE and dif- 
ferent NLTE recipes, the model variation of W in bright regions 
is always steeper than the observations. A higher ybright in the 
models when compared with the observations is more apparent 
for jd < 0.6. 



6. Conclusions 

We have obtained high- spatial-resolution observations of lines 
of neutral oxygen and other species in quiet solar granulation 
at diff'erent positions on the solar disk. These have been used to 
test the 3D photosphere models and line formation. At the so- 
lar disk-centre the predicted distributions of line strengths, shifts 
and shapes are for most lines in excellent agreement with the ob- 
served. Both for the oxygen lines and most of the other lines. At 
diff'erent viewing angles a quantitative study of the line proper- 
ties is more difficult because of a smaller contrast between gran- 
ules and intergranular regions. Nevertheless, fits to the variation 



of line strengths and other quantities over the granulation show 
an encouraging good agreement with the observations. 

For the Oi 777 nm lines we employ the high-spatial- 
resolution observations at disk-centre to constrain the NLTE 
physics. Although not as robust as the centre-to-limb variation 
([ Allende Prieto et aL]|2004|Paper II[ ), this test provides another 
approach to empirically constrain the efficiency of collisions 
with neutral hydrogen, ^h- We find that 5'h = 1 describes most 
observables well, i n agreem ent with the results of [Paper iT Also 
in agreement with Paper II we find that NLTE profiles for = 1 
are narrower than the observations, which agree better with the 
higher Su profiles. Furthermore, it is difficult for any 5'h to de- 
scribe the variation in equivalent width with continuum inten- 
sity in the intergranular regions and granules at the same time. 
These discrepancies may be linked to the crude approximation 
of using scaled classical formulas to estimate the collisions with 
hydrogen. The discrepancies may indicate a more complex de- 
pendence of these collisions with temperature and/or pressure. 
Moreover, Su will also likely vary between transitions. 

Regarding the weak O i 615.81 nm line, its NLTE eff'ects are 
very small and not noticeable in the variation of the line proper- 
ties over the granulation. 

For the spectral lines investigated here, the agreement be- 
tween the 3D model and the observations for the line velocities, 
shapes and strengths is excellent, both for disk-centre and dif- 
ferent viewing angles. It places the 3D model on a solid footing 
regarding oxygen and provides compelling evidence on its suit- 
ability for abundance analysis. 
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